SLAC-PUB-12806 
BAa4R-PUB-07/054 



Study of e+e- AA, AI!°, U^U'^ using Initial State Radiation with BaBAR 

B. Aubert, M. Bona, D. Boutigny, Y. Karyotakis, J. P. Lccs, V. Poircau. X. Prudent, V. Tisscrand, and A. Zghiche 
Laboratoire de Physique des Particules, IN2P3/CNRS et Universite de Savoie, F-74941 Annecy-Le-Vieux, France 

J. Garra Tico and E. Granges 
Universitat de Barcelona, Facultat de Fisica, Departament ECM, E-08028 Barcelona, Spain 

L. Lopez, A. Palano, and M. Pappagallo 
Universita di Bari, Dipartimento di Fisica and INFN, 1-70126 Bari, Italy 

G. Eigen, B. Stugu, and L. Sun 
University of Bergen, Institute of Physics, N-5007 Bergen, Norway 

G. S. Abrams, M. Battaglia, D. N. Brown, J. Button-Shafcr, R. N. Cahn, Y. Groysman, R. G. Jacobsen, 
J. A. Kadyk, L. T. Kerth, Yu. G. Kolomcnsky, G. Kukartsev, D. Lopes Pegna, G. Lynch, L. M. Mir, 
T. J. Orimoto, L L. Osipenkov, M. T. Ronan,* K. Tackmann, T. Tanabe, and W. A. Wenzel 
Lawrence Berkeley National Laboratory and University of California, Berkeley, California 94720, USA 

P. del Amo Sanchez, C. M. Hawkes, and A. T. Watson 
University of Birmingham, Birmingham, B15 2TT, United Kingdom 

H. Koch and T. Schroeder 
Ruhr Universitat Bochum, Institut fiir Experimentalphysik 1, D-44780 Bochum, Germany 

D. Walker 

University of Bristol, Bristol BS8 ITL, United Kingdom 

D. J. Asgeirsson, T. Cuhadar-Donszelmann, B. G. Fulsom, C. Hearty, T. S. Mattison. and J. A. McKenna 
University of British Columbia, Vancouver, British Columbia, Canada V6T IZl 

A. Khan, M. Saleem, and L. Teodorescu 
Brunei University, Uxbridge, Middlesex UB8 3PH, United Kingdom 

V. E. Bhnov, A. D. Bukin, V. P. Druzhinin, V. B. Golubev, A. P. Onuchin, 
S. L Serednyakov, Yu. L Skovpen, E. P. Solodov, and K. Yu. Todyshev 
Budker Institute of Nuclear Physics, Novosibirsk 630090, Russia 

M. Bondioh, S. Curry, L Eschrich, D. Kirkby, A. J. Lankford, P. Lund, M. Mandelkern, E. C. Martin, and D. P. Stoker 

University of California at Irvine, Irvine, California 92697, USA 

S. Abachi and C. Buchanan 
University of California at Los Angeles, Los Angeles, California 90024, USA 

S. D. Foulkes, J. W. Gary, F. Liu, O. Long, B. C. Shcn, G. M. Vitug, and L. Zhang 
University of California at Riverside, Riverside, California 92521, USA 

H. P. Paar, S. Rahatlou, and V. Sharma 
University of California at San Diego, La Jolla, California 92093, USA 

J. W. Berryhill, C. Campagnari, A. Cunha, B. Dahmes, T. M. Hong, D. Kovalskyi, and J. D. Richman 
University of California at Santa Barbara, Santa Barbara, California 93106, USA 



T. W. Beck, A. M. Eisner, C. J. Flacco, C. A. Heusch, J. Kroseberg, W. S. Lockman, 



2 

T. Schalk, B. A. Schumm, A. Seidcn, M. G. Wilson, and L. O. Winstrom 

University of California at Santa Cruz, Institute for Particle Physics, Santa Cruz, California 95064, USA 

E. Chen, C. H. Cheng, F. Fang, D. G. Hitlin, I. Narsky, T. Piatenko, and F. C. Porter 
California Institute of Technology, Pasadena, California 91125, USA 

R. Andreassen, G. Mancinelli, B. T. Meadows, K. Mishra, and M. D. Sokoloff 
University of Cincinnati, Cincinnati, Ohio 45221, USA 

F. Blanc, P. C. Bloom, S. Chen, W. T. Ford, J. F. Hirschauer, A. Kreisel, M. Nagel, 
U. Nauenberg, A. Olivas, J. G. Smith, K. A. Ulmcr, S. R. Wagner, and J. Zhang 
University of Colorado, Boulder, Colorado 80309, USA 

A. M. Gabareen, A. Soffer,t W. H. Toki, R. J. Wilson, and F. Winklmeier 

Colorado State University, Fort Collins, Colorado 80523, USA 

D. D. Altenburg, E. Feltresi, A. Hauke, H. Jasper, J. Merkel, A. Petzold, B. Spaan, and K. Wacker 
Universitdt Dortmund, Institut fiir Physik, D-44221 Dortmund, Cermany 

V. Klose, M. J. Kobel, H. M. Lacker, W. F. Mader, R. Nogowski, 
J. Schubert, K. R. Schubert, R. Schwierz, J. E. Sundermann, and A. Volk 
Technische Universitdt Dresden, Institut fiir Kern- und Teilchenphysik, D-01062 Dresden, Cermany 

D. Bernard, G. R. Bonneaud, E. Latour, V. Lombardo, Ch. Thiebaux, and M. Verderi 
Laboratoire Leprince-Ringuet, CNRS/IN2P3, Ecole Polytechnique, F-91128 Palaiseau, France 

P. J. Clark, W. Gradl, F. Muheim, S. Playfer, A. 1. Robertson, J. E. Watson, and Y. Xie 
University of Edinburgh, Edinburgh EH9 3JZ, United Kingdom 

M. Andreotti, D. Bottom. C. Bozzi, R. Calabrese, A. Cecchi, G. Cibinetto. P. Franchini, 
E. Luppi, M. Ncgrini, A. PctrcUa, L. Piemontese, E. Prcncipe, and V. Santoro 
Universitd di Ferrara, Dipartimento di Fisica and INFN, 1-44100 Ferrara, Italy 

F. AnuUi, R. Baldini-Ferroli, A. Calcaterra, R. dc Sangro. G. Finocchiaro. 
S. Pacctti, P. Patteri, I. M. Pcruzzi,* M. Piccolo, M. Rama, and A. Zallo 
Laboratori Nazionali di Frascati dellTNFN, 1-00044 Frascati, Italy 

A. Buzzo, R. Contri, M. Lo Vetere, M. M. Maori, M. R. Monge, 
S. Passaggio, C. Patrignani, E. Robutti, A. Santroni, and S. Tosi 

Universitd di Genova, Dipartimento di Fisica and INFN, I-I6I46 Cenova, Italy 

K. S. Chaisanguanthum, M. Morii. and J. Wu 

Harvard University, Cambridge, Massachusetts 02138, USA 

R. S. Dubitzky, J. Marks, S. Schenk, and U. Uwer 
Universitdt Heidelberg, PhysikaUsches Institut, Philosophenweg 12, D-69120 Heidelberg, Cermany 

D. J. Bard, P. D. Daunccy, R. L. Flack, J. A. Nash, W. Panduro Vazquez, and M. Tibbetts 
Imperial College London, London, SW7 2AZ, United Kingdom 

P. K. Behera, X. Chai, M. J. Charles, and U. MaUik 
University of Iowa, Iowa City, Iowa 52242, USA 

J. Cochran, H. B. Crawley, L. Dong, V. Eyges, W. T. Meyer, S. Prell, E. I. Rosenberg, and A. E. Rubin 

Iowa State University, Ames, Iowa 5001 1-3160, USA 



Y. Y. Gao, A. V. Gritsan, Z. J. Guo, and C. K. Lae 
Johns Hopkins University, Baltimore, Maryland 21218, USA 



3 



A. G. Dcnig, M. Fritsch, and G. Schott 

Universitdt Karlsruhe, Institut fiir Experimentelle Kernphysik, D-76021 Karlsruhe, Germany 

N. Arnaud, J. Bequilleux, A. D'Orazio, M. Davier, G. Grosdidier, A. Hocker, V. Lepeltier, F. Le Diberder, A. M. Lutz, 
S. Pruvot, S. Rodier, P. Roudeau, M. H. Schune, J. Serrano, V. Sordini, A. Stocchi, W. F. Wang, and G. Wormser 
Laboratoire de I'Accelerateur Lineaire, IN2P3/CNRS et Universite Paris-Sud 11, 
Centre Scientifique d'Orsay, B. P. 34, F-91898 ORSAY Cedex, France 

D. J. Lange and D. M. Wright 

Lawrence Livermore National Laboratory, Livermore, California 94550, USA 

I. Bingham, C. A. Chavez, J. R. Fry, E. Gabathuler, R. Garnet, 
D. E. Hutchcroft, D. J. Payne, K. C. Schofield, and C. Touramanis 
University of Liverpool, Liverpool L69 7ZE, United Kingdom 

A. J. Bcvan, K. A. George, F. Di Lodovico, and R. Sacco 
Queen Mary, University of London, El 4^3, United Kingdom 

G. Cowan, H. U. Flaecher, D. A. Hopkins, S. Paramesvaran, F. Salvatore, and A. C. Wren 
University of London, Royal Holloway and Bedford New College, Egham, Surrey TW20 OEX, United Kingdom 

D. N. Brown and C. L. Davis 
University of Louisville, Louisville, Kentucky 40292, USA 

J. AUison, D. Bailey, N. R. Barlow, R. J. Barlow, Y. M. Chia, C. L. Edgar, G. D. Lafferty, T. J. West, and J. I. Yi 

University of Manchester, Manchester M13 9PL, United Kingdom 

J. Anderson, C. Chen, A. Jawahery, D. A. Roberts, G. Simi, and J. M. Tuggle 
University of Maryland, College Park, Maryland 20742, USA 

G. Blaylock, C. Dallapiccola, S. S. Hertzbach, X. Li, T. B. Moore, E. Salvati, and S. Saremi 
University of Massachusetts, Amherst, Massachusetts 01003, USA 

R. Cowan, D. Dujniic, P. H. Fisher, K. Koeneke, G. ScioUa, 
M. Spitznagel, F. Taylor, R. K. Yamamoto, M. Zhao, and Y. Zheng 
Massachusetts Institute of Technology, Laboratory for Nuclear Science, Cambridge, Massachusetts 02139, USA 

S. E. Mclachhn,* P. M. Patel, and S. H. Robertson 
McGill University, Montreal, Quebec, Canada H3A 2T8 

A. Lazzaro and F. Palombo 
Universitd di Milano, Dipartimento di Fisica and INFN, 1-20133 Milano, Italy 

J. M. Bauer, L. Cremaldi, V. Eschenburg, R. Godang, R. Kroeger, D. A. Sanders, D. J. Summers, and H. W. Zhao 

University of Mississippi, University, Mississippi 38677, USA 

S. Brunet, D. Cote, M. Simard, P. Taras, and F. B. Viaud 
Universite de Montreal, Physique des Particules, Montreal, Quebec, Canada H3C 3J7 

H. Nicholson 

Mount Holyoke College, South Hadley, Massachusetts 01075, USA 

G. De Nardo, F. Fabozzi,^ L. Lista, D. Monorchio, and C. Sciacca 
Universitd di Napoli Federico II, Dipartimento di Scienze Fisiche and INFN, 1-80126, Napoli, Italy 

M. A. Baak, G. Raven, and H. L. Snoek 
NIKHEF, National Institute for Nuclear Physics and High Energy Physics, NL-1009 DB Amsterdam, The Netherlands 



4 



C. P. Jessop, K. J. Knocpfcl, and J. M. LoSecco 

University of Notre Dame, Notre Dame, Indiana 46556, USA 

G. Benelli, L. A. Corwin, K. Honscheid, H. Kagan, R. Kass, J. P. Morris, 
A. M. Rahimi, J. J. Regensburger, S. J. Sekula, and Q. K. Wong 
Ohio State University, Columbus, Ohio 43210, USA 

N. L. Blount, J. Brau, R. Frey, O. Igonkina, J. A. Kolb, M. Lu, 
R. Rahmat, N. B. Sincv, D. Strom, J. Strubc, and E. Torrence 
University of Oregon, Eugene, Oregon 97403, USA 

N. Gagliardi, A. Gaz, M. Margoni, M. Morandin, A. Pompili, 

M. Posocco, M. Rotondo, F. Simonctto, R. Stroili, and C. Voci 

Universita di Padova, Dipartimento di Fisica and INFN, 1-35131 Padova, Italy 

E. Ben-Haim, H. Briand, G. Calderini, J. Chauveau, P. David, L. Del Buono, 
Ch. de la Vaissiere, O. Hamon, Ph. Leruste, J. Malcles, J. Ocariz, A. Perez, and J. Prendki 

Laboratoire de Physique Nucleaire et de Hautes Energies, 
IN2P3/CNRS, Universite Pierre et Marie Curie-Paris6, 
Umversite Denis Diderot-Paris?, F-75252 Pans, France 

L. Gladney 

University of Pennsylvania, Philadelphia, Pennsylvania 19104, USA 

M. Biasini, R. CovarcUi, and E. Manoni 
Universita di Perugia, Dipartimento di Fisica and INFN, 1-06100 Perugia, Italy 

C. Angelini, G. Batignani, S. Bettarini, M. Carpinelli, R. Cenci, A. Cervelli, F. Forti, M. A. Giorgi, 
A. Lusiani, G. Marchiori, M. A. Mazur, M. Morganti, N. Neri, E. Paoloni, G. Rizzo, and J. J. Walsh 
Universita di Pisa, Dipartimento di Fisica, Scuola Normale Superiore and INFN, 1-56127 Pisa, Italy 

J. Biesiada, P. Elmer, Y. P. Lau, C. Lu, J. Olsen, A. J. S. Smith, and A. V. Tclnov 
Princeton University, Princeton, New Jersey 08544, USA 

E. Baracchini, F. Bellini, G. Cavoto, D. del Re, E. Di Marco, R. Faccini, F. Fcrrarotto, F. Fcrroni, M. Gaspero, 
P. D. Jackson, L. Li Gioi, M. A. Mazzoni, S. Morganti, G. Piredda, F. Polci, F. Renga, and C. Voena 
Universita di Roma La Sapienza, Dipartimento di Fisica and INFN, 1-00185 Roma, Italy 

M. Ebcrt, T. Hartmann, H. Schroder, and R. Waldi 
Universitdt Rostock, D-18051 Rostock, Germany 

T. Adye, G. Castelh, B. Franek, E. O. Olaiya, W. Roethel, and F. F. Wilson 
Rutherford Appleton Laboratory, Chilton, Didcot, Oxon, 0X11 OQX, United Kingdom 

S. Emery, M. Escalier, A. Gaidot, S. F. Ganzhur, G. Hamel de Monchenault, 
W. Kozanecki, G. Vasseur, Ch. Yeche, and M. Zito 
DSM/Dapma, CEA/Saclay, F-91191 Gif-sur-Yvette, France 

X. R. Chen, H. Liu, W. Park, M. V. Purohit, R. M. White, and J. R. Wilson 
University of South Carolina, Columbia, South Carolina 29208, USA 

M. T. Allen, D. Aston, R. Bartoldus, P. Bechtle, R. Glaus, J. P. Coleman, M. R. Convery, J. C. Dingfelder, 
J. Dorfan, G. P. Dubois-Felsmann, W. Dunwoodie, R. C. Field, T. Glanzman, S. J. Gowdy, M. T. Graham, 
P. Grenier, C. Hast, W. R. Innes, J. Kaminski, M. H. Kelsey, H. Kim, P. Kim, M. L. Kocian, D. W. G. S. Leith, 
S. Li, S. Luitz, V. Luth, H. L. Lynch, D. B. MacFarlanc, H. Marsiske, R. Messner, D. R. MuUer, C. P. O'Grady, 
L Ofte, A. Perazzo, M. Perl, T. Pulham, B. N. Ratchff, A. Roodman, A. A. Salnikov, R. H. Schindler, 
J. Schwiening, A. Snyder, D. Su, M. K. Sullivan, K. Suzuki, S. K. Swain, J. M. Thompson, J. Va'vra, A. P. Wagner, 



5 



M. Weaver, W. J. Wisnicwski, M. Wittgen, D. H. Wright, A. K. Yarritu, K. Yi, C. C. Young, and V. Ziegler 

Stanford Linear Accelerator Center, Stanford, California 94309, USA 

P. R. Burchat, A. J. Edwards, S. A. Majewski, T. S. Miyashita, B. A. Petersen, and L. Wilden 

Stanford University, Stanford, California 94305-4060, USA 

S. Ahmed, M. S. Alam, R. Bula, J. A. Ernst, V. Jain, B. Pan, M. A. Saeed, F. R. Wappler, and S. B. Zain 

State University of New York, Albany, New York 12222, USA 

M. Krishnamurthy and S. M. Spanier 

University of Tennessee, Knoxville, Tennessee 37996, USA 

R. Eckmann, J. L. Ritchie, A. M. Ruland, C. J. Schilhng, and R. F. Schwitters 
University of Texas at Austin, Austin, Texas 78712, USA 

J. M. Izcn, X. C. Lou, and S. Ye 

University of Texas at Dallas, Richardson, Texas 75083, USA 

F. Bianchi, F. Gallo, D. Gamba, and M. PeUiccioni 
Universitd di Torino, Dipartimento di Fisica Sperimentale and INFN, 1-10125 Torino, Italy 

M. Bombcn, L. Bosisio, C. Cartaro, F. Cossutti, G. DcUa Ricca, L. Lanccri, and L. Vitalc 
Universitd di Trieste, Dipartimento di Fisica and INFN, 1-34127 Trieste, Italy 

V. Azzohni, N. Lopez-March, F. Martinez- Vidal,^ D. A. Milanes, and A. Oyanguren 
IFIC, Universitat de Valencia-CSIC, E-46071 Valencia, Spam 

J. Albert, Sw. Banerjee, B. Bhuyan, K. Hamano, R. Kowalewski, L M. Nugent, J. M. Roney, and R. J. Sobie 

University of Victoria, Victoria, British Columbia, Canada V8W 3P6 

P. F. Harrison, J. Ihc, T. E. Latham, and G. B. Mohanty 
Department of Physics, University of Warwick, Coventry CV4 7AL, United Kingdom 

H. R. Band, X. Chen, S. Dasu, K. T. Flood, J. J. Hollar, P. E. Kutter, Y. Pan, M. Pierini, R. Prcpost, and S. L. Wu 

University of Wisconsin, Madison, Wisconsin 53706, USA 

H. Neal 

Yale University, New Haven, Connecticut 06511, USA 

We study the e'^e~ AA'y, AS^'y, S^S'^'j processes using 230 fb^^ of integrated luminosity 
collected by the BABAR detector at e^e~ center-of-mass energy of 10.58 GeV. From the analysis of 
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In this paper we continue the experimental study of 
baryon time- like electromagnetic form factors. In our 
previous work [1] we have measured the energy depen- 
dence of the cross section for e"^e~ pp and of the 
proton form factor using the initial state radiation (ISR) 
technique. Here we use this technique to study the pro- 
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FIG. 

e+e- 



yasR) 



/ 



The Feynman diagram describing the ISR process 
/7, where / is a hadronic system. 



> AA, AS". The Born cross sec- 

tion for the ISR process e+e~ — > f + j (Fig.l). where / is 
a hadronic system, integrated over the hadron momenta, 
is given by 



dcT, 



dm d cos 9* 



2m , 
W{x, 



(1) 



where is the e~^e center-of-mass energy (cm.), m 
is the invariant mass of the hadronic system, af{m) is 



the cross section for 



/ reaction. 



1 — rr? I and and 6'* arc the ISR photon energy and 
polar angle, respectively, in the e^e^ cm. frame. ^ The 
function [2] 



baryon momentum in the dibaryon rest frame and the 
momentum of the BB system in the e~ cm. frame. 
This distribution can be expressed as the sum of the 
terms proportional |Gj\/p and IG^p. The Qb depen- 
dencies of the Ge and Gu terms are close to sin^ Qb and 
1 + cos^ Ob angular distributions for electric and mag- 
netic form factors in the e+e^ BB process. The full 
differential cross section for e^e~ BB^ [3] is given in 
the Appendix. 

A nonzero relative phase between the electric and mag- 
netic form factors manifests itself in polarization of the 
outgoing baryons. In the e^e~ —> BB reaction this po- 
larization is perpendicular to the production plane [4]. 
For the ISR process e+e~ —>■ BBj the polarization ob- 
servables are analyzed in Refs. [3, 5]. The expression for 
the baryon polarization as a function of Ge, Gm, and 
momenta of the initial electron, ISR photon, and final 
baryon [3] is given in the Appendix. In the case of the 
AA final state the A — *■ pTT~ decay can be used to mea- 
sure the A polarization and hence the phase between the 
form factors. 

Experimental information on the e^e^ AA, S'^S", 
AS'^ reactions is very scarce. The e+e" — > AA cross 
section is measured as 100^35 pb at 2.386 GeV, and at 
the same energy upper limits for e+e~ — > S'^S" (< 120 
pb) and e+e" ^ AS'^ (< 75 pb) cross sections have been 
obtained [6]. No other experimental results exist. 



a 



2-2x + x'^ 



nx \ sin^ 61* 



(2) 



II. THE BABAR DETECTOR AND DATA 
SAMPLES 



describes the probability of ISR photon emission for 9* ^ 
nifilsfs, where a is the fine structure constant and me 
is the electron mass. The cross section for the process 
e+e" — > BB, where is a spin-1/2 baryon, depends on 
magnetic (Gm) and electric [Ge) form factors as follows: 



' BB 



(m) 



|GM(m)p + ^|Gs(m)p 



(3) 



where /? ~ -^/l — 477i^ / m^ and r = m^/4m^; at thresh- 
old, Ge = Gm- The cross section determines the linear 
combination of the squared form factors 



1^.^,2 ^ 2r|GM(m)p + |G^(m)p 



(4) 



and we define |F(m)| to be the effective form factor [1]. 

The modulus of the ratio of electric and magnetic form 
factors can be determined from the analysis of the dis- 
tribution of cos 9 B , where is the angle between the 



^ Throughout this paper the use of charge conjugate modes is im- 
plied. 

^ Throughout this paper the asterisk denotes quantities in the 
e+e~ cm. frame. 



We analyse a data sample corresponding to an inte- 
grated luminosity of 230 fb~^ recorded with the BABAR 
detector [7] at the PEP-II asymmetric-energy storage 
rings. At PEP-II, 9-GeV electrons collide with 3.1-GeV 
positrons at a center-of-mass energy of 10.58 GcV (the 
T(4S) resonance). Additional data (~ 10%) recorded at 
10.54 GeV are included in the present analysis. 

Charged-particle tracking is provided by a five-layer 
silicon vertex tracker (SVT) and a 40-layer drift chamber 
(DCH), operating in a 1.5-T axial magnetic field. The 
transverse momentum resolution is 0.47% at 1 GeV/c. 
Energies of photons and electrons are measured with a 
CsI(Tl) electromagnetic calorimeter (EMC) with a reso- 
lution of 3% at 1 GeV. Charged-particle identification is 
provided by specific ionization {dE/dx) measurements in 
the SVT and DCH, and by an internally reflecting ring- 
imaging Cherenkov detector (DIRC). Muons are identi- 
fied in the solenoid's instrumented flux return, which con- 
sists of iron plates interleaved with resistive plate cham- 
bers. 

Signal ISR processes are simulated with the Monte 
Carlo (MC) event generator Phokhara [8, 9]. Because 
the polar-angle distribution of the ISR photon is peaked 
near 0° and 180°, the MC events are generated with a 
restriction on the photon polar angle: 20° < 9* < 160°. 
The Phokhara event generator includes next-to-lcading- 
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order radiative corrections to the Born cross section. In 
particular, it generates an extra soft photon emitted from 
the initial state. To restrict the maximum energy of the 
extra photon we require that the invariant mass of the 
dibaryon system and the ISR photon satisfies M^g^ > 8 
GeV/c^. The generated events are subjected to detailed 
detector simulation based on GEANT4 [10], and are re- 
constructed with the software chain used for the experi- 
mental data. Variations in the detector and in the beam 
background conditions arc taken into account. For the 
full simulation we use the differential cross section for 
the e~^e~ BBj process with Ge — Gm- In order 
to study angular distributions and model dependence of 
detection efficiency we produce two large samples of sim- 
ulated events at the generator level, one with = 
and the other with Gm = 0, and reweight the events 
from the full simulation sample according to the desired 
\Ge/Gm\ ratio. 

Background from e^e^ — s- qq, where q represents a 
M, d, s or c quark, is simulated with the JETSET [11] 
event generator. JETSET also generates ISR events with 
hadron invariant mass above 2 GeV/c^ and therefore can 
be used to study ISR background with baryons in the 
final state. The most important background processes 
e+e" ^ BBtt^j, e+e" BBtt°, and e+e" ApK+j 
are simulated separately. Three-body phase space and 
the Bonneau-Martin formula [2] are used to generate the 
angular and energy distributions for the final hadrons 
and ISR photon, respectively. For these processes extra 
soft-photon radiation from the initial state is generated 
using the structure fimction method [12]. 

III. THE REACTION e+e" AA-y 
A. Event selection 

The initial selection of events requires the presence of 
a high energy photon and at least one A and one A can- 
didate. The hard photon must have energy in the cm. 
frame E* > 3 GeV. The A — > pTr~ decay mode with the 
branching fraction of (63.9 ± 0.5)% [13] is used to iden- 
tify A candidates. Two oppositely-charged tracks are as- 
signed the proton and pion mass hypotheses and fitted to 
a common vertex. Any combination with invariant mass 
in the range 1.104-1.128 GeV/c^ (the nominal A mass is 
1.115683(6) GeV/c^ [13]), laboratory momentum greater 
than 0.5 GeV/c, and fit probability greater than 0.001 
is considered a A candidate. The candidate is then re- 
fitted with a A mass constraint to improve the precision 
of the A momentum measurement. To suppress combi- 
natorial background we require that at least one of the 
proton candidates be identified as a proton according to 
the specific ionization {dE/dx) measured in the SVT and 
DCH, and the Cherenkov angle measured in the DIRC. 

For events passing the preliminary selection, we per- 
form a kinematic fit that imposes energy and momentum 
conservation at the production vertex to the A and A can- 



didates and the photon with highest E* . For events with 
more than one A (A) candidate we consider all possible 
ylyl combinations, and the one giving the lowest for 
the kinematic fit is retained. The MC simulation does not 
accurately reproduce the shape of the resolution function 
for the photon energy. This leads to a difference in the 
distributions resulting from the kinematic fits to data 
and simulated events. To reduce this difference, only the 
measured direction of the ISR photon is used in the fit; 
its energy is a fit parameter. The distributions for 
the kinematic fit (x'aa) ^o data events and to simulated 
e+e^ AAj events arc shown in Fig. 2. We select the 
events with x'aa < 20 for further analysis. The control 
region 20 < x'aa used for background estimation 

and subtraction. 

Possible sources of background for the process under 
study are those with only one A in the final state, such 
as e+e" ApK~^-f. Such events contain a charged kaon 
instead of one of the pion candidates. To suppress this 
background we require that no charged pion candidate 
be identified as a kaon. This requirement rejects 70% of 
the background from e+e" — > ApK^j and only ^2% of 
signal events. 

The scatter plot of the invariant mass of the A candi- 
date versus the invariant mass of the A candidate for the 
387 data events passing all the selection criteria is shown 
in Fig. 3(a) and that for simulated e+e~ AAj events 
is shown in Fig. 3(b). The AA invariant mass spectrum 
for data events is shown in Fig. 4. About half of the 
events have invariant mass below 3 GeV/c^. Signals due 
to J/ip AA and ip{2S) AA decays are also clearly 
seen. 



B. Background subtraction 

Processes of three kinds potentially contribute back- 
ground to the e'^e" — > AAj data sample, namely those 
with zero, one and two yl's in the final state. 

The composition of the one-yl background is studied 
using JETSET simulation. This background is domi- 
nated by e+e^ ApK^^ events. Other processes also 
contain a charged kaon in the final state. The level of the 
one-Zl background can be estimated from the fraction of 
data events rejected by the requirement that no 7r+ can- 
didate be identified as a K+ . For M^j < 3 GeV/c^ this 
fraction is 3/224, and we estimate that the one-yl back- 
ground does not exceed 1.6 events at 90% confidence level 
(CL). 

A more precise estimation (but based on JETSET pre- 
diction for the composition of onc-yl events) of this back- 
ground is obtained using a special selection of e^e~ 
ApK~^j events. We select events with at least 4 charged 
tracks and a photon with E* > 3 GeV. Two tracks, 
one of which is identified as a proton, must be com- 
bined to form a A candidate, and the other two must 
originate from the e~^e~ interaction point and be iden- 
tified as an antiproton and a positively-charged kaon. 
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FIG. 2: The Xaa distributions for data (a) and e+e — > AA'j simulation(b). 
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FIG. 3: Scatter plots of the invariant mass of the A candidate versus the invariant FIG. 4: The AA mass spectrum for 
mass of the A candidate for data (a) and e"'"e~ — > AA'j simulation (b). events satisfying the AAy selection cri- 



teria. 



For these events we perform the kinematic fit to the 
e^e~ — > ApK^j hypothesis and require < 20. The 
background for e"*"e~ ApK^^ is estimated from the re- 
gion 20 < < 40. The total number of selected ApK~^j 
events is found to be 568 ± 30. Using the ratio of detec- 
tion efficiencies for AA'j and ApK~^j selections obtained 
from simulation, (0.12±0.07)%, we calculate the number 
of ApK~^j events satisfying the AA^y selection criteria to 
be 0.7 ± 0.4. Taking the ratio of ApK'^j to all one-yl 
events (0.7) from JETSET simulation we estimate the 
total number of one-vl background events to be 1.0 ±0.6. 
The e'^e^ ApK'^^ simulation is reweighted to repro- 
duce the shape of the experimental M^pK distribution. 
The reweighted events are then used to find the distribu- 
tion of M^j for events with only one real A. We estimate 
0.8 ± 0.5 background events to have Af^j < 3 GeV/c^. 

The background processes with no real yl's are 
ISR processes with four charged particles in the final 
state: e+e~ 27r+27r~7, K'^K~tt^-k~^^ K^Kstt~j, 
pp7r"''7r~7, 27r"'"27r~7r°, etc. The background from these 
processes can be estimated from an analysis of the two- 



dimensional distribution of the A and A candidate mass 
values. The 6x6 two-dimensional histogram correspond- 
ing to the plot in Fig. 3a is fitted by the following func- 
tion: 

= N2S,Sj + NqBq,Boj + Ni{S,Bij + S,Bu)/2. (5) 

where the six mass intervals are those shown in Fig. 5. 
Here A''2, A^i, and Nq represent the numbers of events 
with two, one, and zero yl's, respectively; A^2 and Nq are 
free fit parameters, and Ni is fixed at the value deter- 
mined above (0.8 ± 0.5 events for M^^ < 3 GcV/c^); S^ 
is the probability for a yl to have reconstructed mass in 
the ith mass bin, while Boi and Bn are the probabilities 
for a false A candidate from background with zero and 
one real A, respectively. Since the one-yl background is 
small and its presence leads to only small changes in the 
fitted N2 and Nq values, we use a uniform distribution for 
Bi(m), i.e. all Bu = 1/6. The Bo{m) are parametrized 
by the linear function Bqi = 1/6 + (i — 3. 5) A, and five 
of the Si and A arc free fit parameters. For 221 events 
with Af^j < 3 GeV/c2 the fit yields N2 = 216 ± 15 
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FIG. 5: The values of Si (see text) obtained from fits to 
data (points with error bars) and from fits to e'^e~ — > AA'y 
simulated events (histogram). 
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FIG. 6: The distribution of data events satisfying the AA'y 
selection criteria over chosen mass intervals. The shaded his- 
togram shows fitted background. 



and No = itj- The fitted values of the Si are m good 
agreement with the values expected from e+e^ AAj 
simulation (Fig. 5). In particular, S3 + S4 = 0.950±0.014 
for data and 0.953 ± 0.003 for simulation. 

The sources of two-yl background are processes with 
extra neutral particle(s) in the final state: e+e~ AAtt'^, 
e+e" AA'y^, e+e" AAtt'^'j, etc. A significant frac- 
tion of e+e^ —>■ AAtt'^ events with an undetected low- 
energy photon or with merged photons from tt*^ decay 
are reconstructed under the AA'y hypothesis with a low 
value of XaA' ''^^'^ ^^ot t)e separated from the process 
under study. This background is studied by selecting a 
special subsample of events containing a A and a A can- 
didate and at least two photons, one with energy greater 
than 0.1 GeV and the other with cm. energy above 3 
GeV. The two-photon invariant mass is required to be 
in the range 0.07 to 0.2 GeV/c^ A kinematic fit to the 
e~^e~ AA'yj hypothesis is then performed. Require- 



ments on the (x^ < 20) and the two-photon invariant 
mass (0.11 < M^^ < 0.16 GeV/c^) are imposed to de- 
fine AAtt'^ candidates. No data events satisfy these crite- 
ria, and the expected background from e'^e~ —f AA'y is 
0.8 ± 0.3 events. The corresponding 90% CL upper limit 
on the number of selected AAtt'^ candidates is 1.6 events. 
Using the ratio of detection efficiencies for AAtt^ and AA'y 
selections (0.28±0.02) we find that the AAtt'^ background 
in AA'y sample does not exceed 6 events. This upper limit 
is used as a measure of the systematic uncertainty due to 
AA'K^ background. We assume that the dibaryon mass 
spectrum in the e'^e~ — > AAn'^ process is similar to that 
for the pp-K^ final state [1]. In particular, about 70% of 
ylyl7r° events are located in the AA mass region below 3 
GeV/c^. For this mass range this background does not 
exceed 2% of the selected AA'y candidates. 

The two-yl background other than from e+e~ — > AAtt'^ 
has the Xaa distribution very different from that for the 
process under study. Tabic I shows the ratio of num- 
bers of selected ylyl7 candidates with 20 < Xaa < 40 
and Xaa < 20 for signal and background processes. The 
ratios are obtained from simulation. The column de- 
noted "JETSET" shows the result of JETSET simula- 
tion for background events containing two yl's in the fi- 
nal state. From the number of selected two-vl events 
in the signal and control Xaa regions, N2{x^ < 20) and 
A^2(20 < < 40), the numbers of signal and background 
events with Xaa < 20 can be calculated as: 

^ A^2(x' < 20) - jV2(20 < < 40)/Pbkg 

1 — Psig/Pbkg 

N2b = iV2(x' < 20)-7V2., (6) 

where (3bkg is the ratio of fractions of events in the con- 
trol and signal x^ regions averaged over all processes 
contributing into two-Zl background. For this coefficient 
we use f3bkg = 0.9 ± 0.3 which is close to the value ob- 
tained from the JETSET simulation, with an uncertainty 
covering the /?; variations for different background pro- 
cesses. For the ratio for the signal process Psig, we use 
the value obtained from simulation /3sig = 0.073 ± 0.010. 
The quoted error takes into account MC statistics, the 
data-MC simulation difference in x^ distribution, and the 
f3sig variation as a function of AA mass. The difference 
between data and simulated x^ distributions was stud- 
ied in Ref. [1] using the process e^e~ — > jX^ jjT^. The 
resulting values of N2S and A^26 for AA masses below 3 
GeV/c^ are listed in Table II. The total background in 
the signal x^ region from the processes with zero, one, 
and two yl's in final state is about 10%. The last row 
of the table shows the JETSET prediction for signal and 
background events in the signal x^ region. The simula- 
tion overestimates the signal yield, but can be used for 
qualitative estimation of background level. 

The procedure for background estimation and subtrac- 
tion described above is applied in each of the twelve AA 
mass intervals indicated in Table IV. Due to restricted 
statistics we fit the two-dimensional histogram of M^^- 
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TABLE L The Pi values obtained from simulation for signal and background processes, where /3i is the ratio of the number of 
selected AA-y candidates with 20 < x\a < 40 to that with xaa < 20. 



AA-y yir"7 S°i:°j AS°-k" jetset 

/3i 0.073 ±0.005 0.83 ±0.07 1.1 ±0.2 0.81 ± 0.09 0.86 ± 0.06 



TABLE II: A'^ is the number of selected AA'y candidates with My^j < 3 GeV/c^, is the number of signal events, Nq, 
Ni, N2b indicate the number of background events with zero, one, and two yl's in the final state, respectively, and N^j^^o is 
background from e+e~ — > AAiv''^. 





N 


N2s 


No 


Ai 


N2t 


Naa.o 


Xaa < 20 


221 


204 ± 19 




0.8 ±0.5 


12± 10 


< 4 


20 < xIa < 40 


35 


15 ±3 




0.6 ±0.4 


11 ±8 


< 1 


xIa < 20 (JETSET) 


522 


500 ± 17 


2.5 ± 1.2 


0.6 ±0.6 


18 ±3 


1.2 ±0.9 



vs Mp^+ using 3x3 bins, and fix the Si (sec Eq.(5)) at 
the values obtained from MC simulation. The histograms 
for signal and control regions arc fitted simultaneously. 
The free fit parameters are Nq and A for the two re- 
gions, N2s, and N2b- Fig. 6 shows the distribution of se- 
lected events over the chosen mass intervals. The shaded 
histogram shows the background contribution obtained 
from the fit. The resulting numbers of signal events are 
listed in Tabic IV, where the quoted errors include the 
statistical errors and errors due to uncertainties in the 
Psig, Pbkg and Si Coefficients. These coefficients are var- 
ied within their uncertainties during fitting. For the mass 
ranges 3.2 < M^^ < 3.6 GeV/c^ and 3.8 < M^j < 5.0 
GcV/c^ where we do not see evidence for a signal above 
background, 90% CL upper limits on the number of sig- 
nal events are listed. The mass regions near the J/ -0 and 
ip{2S) will be considered separately in Sec. IIIF. 



C. Angular distributions for e+e 



AA'y 



The modulus of the ratio of the electric and magnetic 
form factors can be extracted from an analysis of the 
distribution of cos 9a, where 9a is the angle between the 
A momentum in the AA rest frame and the momentum of 
the A A system in the e~^e~ cm. frame. This distribution 
is given by 



diV 



d cos 9 a 



A[Hm{cos9a,M^j) 



G 



M 



He{cos9a,Maa)\- 



(7) 



The functions Hm{cos9a, M^a) and He{cos9a, M^a) do 
not have an analytic form, and so are calculated using 
MC simulation. To do this two samples of e+e^ AA'y 
events were generated, one with Ge = and the other 
with Gm = 0, using generator level simulation. The an- 
gular dependencies of the resulting functions do not differ 
significantly from the (1 ± cos^ 9a) and sin^ 9a functions 



corresponding to the magnetic and electric form factors 
in the case of e+e^ AA. 

The observed angular distributions arc fitted in two 
mass intervals: from AA threshold to 2.4 GeV/c^ and 
from 2.4 GeV/c^ to 2.8 GeV/c^. For each mass and an- 
gular interval, the background is subtracted by means of 
the procedure described in the previous section. The an- 
gular distributions obtained arc shown in Fig. 7. The 
distributions are fitted using the expression on the right- 
hand side of Eq. (7) with two free parameters A and 
\Ge/Gm\- The functions i?M and TJ^; are replaced by the 
histograms, obtained from MG simulation with the AA 
selection criteria applied. To take into account the effect 
of these criteria (Fig. 8), the simulated events produced 
assuming Ge = Gm are reweighted according to the 
cos^yi distributions obtained at generator level. These 
weight functions also take account of the difference in AA 
mass dependence between data and MC simulation. The 
histograms fitted to the angular distributions are shown 
in Fig. 7. The following values of the \Ge/Gm\ ratio are 
obtained: 



\Ge/G 



1.73: 



-0.99 
-0.57 



\Ge/G 



0.71 



+0.66 
-0.71 



(2.23-2.40 GcV)/c2 



(2.40-2.80 GcV)/c2. 



The quoted errors include both statistical and system- 
atic uncertainties. The net systematic uncertainty does 
not exceed 15% of the statistical error and includes the 
uncertainties due to background subtraction, limited MC 
statistics, and the mass dependence of the \Ge/Gm\ ra- 
tio. 

We also measure the angular distribution for J/rp 
AA decay, for which the shape is usually described by 
the form (1 ± acos^6'). The world average value of 
a = 0.65 ± 0.10 [14-16]. The distribution for J/ip AA 
decay in the present experiment is shown in Fig. 9. To 
remove background, this distribution was obtained as the 
difference between the histogram for the signal mass re- 
gion (3.05-3.15 GcV/c^) and that for the mass sidebands 
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FIG. 7: The cos6'a distribution for the mass regions 2.23-2.40 GeV/c^ (a), and 2.40-2.80 GeV/c^ (b). The points with 
error bars represent the data after background subtraction. The histograms are fit results: the dashed histogram shows the 
contributions corresponding to the magnetic form factor; the dotted histogram shows the contribution from the electric form 
factor, and the solid histogram is the sum of these two. 
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FIG. 8: The cos 9a dependence of the detection efficiency for 
simulated events with M^j < 2.8 GeV/c^. 
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FIG. 9: The cos Oa distributions for data (points with error 
bars) and simulation (histogram) corresponding to the reac- 
tion e^e~ J/i^'y AA'y. 



FIG. 10: The ylyl mass dependence of detection efficiency 
obtained from MC simulation. 

(3.00-3.05 and 3.15-3.20 GeV/c^). The data distribution 
is in good agreement with that obtained from simulation 
with a = 0.65. 

Our results on the \Ge/Gm\ ratio arc consistent both 
with \Ge/Gm\ = 1, valid at the AA threshold, and with 
our results for the reaction e^e~ — > pp for which this 
ratio was found to be greater than unity near thresh- 
old [1]. The strong dependence of the \Ge/Gm\ ratio on 
the dibaryon mass near threshold is expected due to the 
baryon-antibaryon final state interaction [17, 18]. 

D. Mass dependence of the detection efficiency 

To first approximation the detection efficiency is de- 
termined from MC simulation as the ratio of true AA 
mass distributions computed after and before applying 
the selection criteria. Since the e+e" — > AAj differential 
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FIG. 11: The distribution of A flight length for data (points 
with error bars) and e'^e^ AA'y simulation (histogram). 



TABLE III: The values of the various efficiency corrections 
Si for the process e'^e~ AA'y. 



effect 


-5., (%) 


Xaa < 20 


-2.0 ±2.0 


no identified K 


-1.0 ± 1.3 


track reconstruction 


-1.0 ±3.8 


p nuclear interaction 


±1.0 ±0.4 


PID 


±0.6 ±0.6 


photon inefficiency 


-1.3 ±0.3 


photon conversion 


±0.4 ±0.2 


trigger 


-0.6 ± 0.5 for M^j < 2.4 GeV/c^ 



total -3.9 ± 4.6 for M^j < 2.4 GeV/c^ 

-3.3 ± 4.6 for M^j > 2.4 GeV/c^ 



cross section depends on two form factors the detection 
efficiency cannot be determined in a model-independent 
way. We use a model in which the \Ge/Gm\ ratio is set 
to the values obtained from the fits to the experimental 
angular distributions for M^j < 2.8 GeV/c^, and then 
set \Ge/Gm\ — 1 for higher masses. The detection ef- 
ficiency obtained in this way is shown in Fig. 10. This 
efficiency includes the branching fraction for A —f p7r" 
decay, which is (63.9±0.5)% [13]. For M^j < 2.8 GeV/c^ 
the variation of the \Ge/Gm\ ratio within its experimen- 
tal uncertainties leads to a 2.5% model uncertainty. For 
higher masses, the model uncertainty is taken as half the 
difference between the detection efSciencies correspond- 
ing Ge = and Gm = 0; this yields a 5% uncertainty. 

The efficiency determined from MC simulation (smc) 
must be corrected to account for data-MC simulation dif- 
ferences in detector response: 

e^SMcYla + S,), (8) 

i 

where the Si's correct for the several effects discussed 
below, and summarized in Table III. 

The efficiency correction for the requirement was 



studied in Ref. [1] for e+e 7 and in Ref. [19] 

for e'^e^ 27r+27r^7. The corrections were found to be 
-(1.0 ± 1.3)% and -(3 ±2)%, respectively. For e+e" ^ 
AA'y we double the correction for fi'^ and assign a 
systematic uncertainty equal to the correction. 

The effect of requiring no identified K is studied using 
e'^e~ —> J/ipj — ^ AA'y events. The number of J/?/; events 
is determined using the sideband subtraction method. 
The event losses when requiring no identified K are found 
to be (2.1 ± 1.2)% in data and (1.1 ± 0.4)% in MC sim- 
ulation. The difference of these numbers is taken as the 
efficiency correction. 

Another source of data-MC simulation difference is 
track loss. The correction due to the difference in track 
reconstruction is estimated to be -0.25% per track with 
systematic uncertainty 0.7% for each proton and 1.2% 
for each pion, which has a softer momentum spectrum. 
Specifically, for the antiproton track only, an extra sys- 
tematic error originates from imperfect simulation of nu- 
clear interactions of antiprotons in the detector material. 
This effect was studied in [1], and the corresponding effi- 
ciency correction is found to be (1.0 ± 0.4)%. All correc- 
tions for track reconstruction described above were ob- 
tained for tracks originating from the e^e~ interaction 
point. To estimate possible data-MC simulation differ- 
ence due to A flight path we compare the distributions 
of reconstructed A flight length (Fig. 11). The data and 
simulated distributions are in good agreement, and so 
there is no need to introduce an extra efficiency correc- 
tion for this effect. 

The data-MC simulation difference for proton identi- 
fication is calculated using the p/p identification proba- 
bilities for data and simulation obtained in Ref. [1] for 

A correction must be also applied to the photon detec- 
tion efficiency. There are two main sources for this cor- 
rection: data-MC simulation difference in the probability 
of photon conversion in the detector material before the 
DCH, and the effect of dead calorimeter channels. Both 
effects were studied in Ref. [1] using e'^e^ //+/i^7 and 
e+e^ — > 77 events. 

The quality of the simulation of the trigger efficiency 
was also studied. The overlap of the samples of events 
passing different trigger criteria, and the independence of 
these triggers, were used to measure the trigger efficiency. 
A small difference (—(0.6 ±0.5)%) in trigger efficiency 
between data and MC simulation was observed for AA 
masses below 2.4 GeV/c^. 

The total efficiency correction is —(3.9 ± 4.6)% for 
M^j < 2.4 GeV/c2 and -(3.3 ± 4.6)% for M^j > 2.4 
GeV/c^. The corrected detection efficiencies are listed 
in Table IV. The uncertainty in detection efficiency in- 
cludes a simulation statistical error, a model uncertainty, 
the error on the A p7r~ branching fraction, and the 
uncertainty of the efficiency correction. 
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TABLE IV: The AA invariant mass interval (My^j), net number of signal events (Ns), detection efficiency (e), ISR luminosity 
(L), measured cross section (a), and effective form factor (F) for e+e^ — > AA. The quoted errors on a are statistical and 
systematic, respectively. For the form factor, the total error is listed. 







e 


L 




\F\ 


(GeV/c2) 






(Pb-^) 


(p'b) 





2.23-2.27 


99 0+6.7 


0.055 ± 0.006 


1.98 


2041^^ ± 22 


U.^OO_0 Q44 


2.27-2.30 


24.31^° 


0.057 ± 0.005 


2.10 


202tl° ± 18 


n 107+0.025 


2.30-2.35 


32.6+t'l 


0.068 ± 0.005 


3.06 


155125 ± 12 


n 1 r.A+0.014, 


2.35-2.40 


35.61^:^3 


0.064 ± 0.005 


3.14 


1761^1 ± 15 


U.10Z_Q QJ^g 


2.40-2.45 


19.2tli 


0.066 ± 0.006 


3.22 


SOlao ± 8 


+0.017 
U.iUO_o.o20 


2.45-2.50 


21.4lt« 


0.062 ± 0.006 


3.30 


104l^t ± 10 


Q -,10+0.013 


2.50-2.60 


22.3+li 


0.070 ± 0.005 


6.85 


46l " ± 4 


u.urz_o 008 


2.60-2.70 




0.064 ± 0.005 


7.18 


251^3 ± 2 


U.UOZ_Q 016 


2.70-2.80 


4 7+4-0 

' -3.6 


0.063 ± 0.006 


7.52 


lOlg ± 1 


f^ f,nO + 0.012 


2.80-3.00 


9 9+3.4 
^•^-3.7 


0.065 ± 0.006 


16.09 


2.ll35 ±0.2 


016+'^ °°^ 


3.20-3.60 


< 4.6 


0.055 ± 0.005 


39.88 


<'2.1 


< 0.017 


3.80-5.00 


< 3.9 


0.066 ± 0.005 


180.38 


< 0.3 


< 0.009 



E. Cross section and form factor 

The cross section for e+e^ AA is calculated from 
the AA mass spectrum using the expression 



a(m) 



(diV/dm)e. 
eRdL/dr. 



(9) 



where {dN/d'm)corr is the mass spectrum corrected for 
resolution effects, dL/dm is the so-called ISR differential 
luminosity, e is the detection efficiency as a function of 
mass, and i? is a radiative correction factor accounting 
for the Born mass spectrum distortion due to emission of 
extra photons by the initial electron and positron. The 
ISR luminosity is calculated using the total integrated 
luminosity L and the probability density function for ISR 
photon emission (Eq. (2)): 

dL OL [ , 9,,1 + C 2to , , 

— = — 2 - 2x + In - x^C L. 10 

dm Tra; \ 1 — C } s 

Here C — cos 0g , and 0g determines the range of polar 
angles of the ISR photon in the e+e" cm. frame: 0g < 
0* < 180° - In our case Qq is equal to 20°, since we 
determine efficieney using simulation with 20° < 0* < 
160°. The values of ISR luminosity integrated over the 
corresponding mass interval are listed in Table IV. 

The radiative correction factor i? is determined us- 
ing Monte Carlo simulation (at the generator level, with 
no detector simulation). The ylyl mass spectrum is 
generated using only the pure Born amplitude for the 
e+ e~ — > AA^ process, and then using a model with next- 
to-leading-order radiative corrections included. The ra- 
diative correction factor, evaluated as the ratio of the 
second spectrum to the first, is found to be practically in- 
dependent of mass, with an average value equal to 1.0035 
for masses below 3 GeV/c^. It should be noted that the 



value of i? depends on the criterion applied to the invari- 
ant mass of the AA^ system. The value of R obtained 
in our case corresponds to the requirement M^j^^ > 8 
GeV/c^ used in our simulation. The theoretical uncer- 
tainty in the radiative correction calculation is estimated 
to be less than 1% [8]. The calculated radiative correc- 
tion factor does not take into account vacuum polariza- 
tion, and the contribution of the latter is included in the 
measured cross section. 



The dependence of the mass resolution on the AA in- 
variant mass is shown in Fig. 12. The mass resolution 
is calculated in simulation as the RMS deviation of the 
A'/^j — M''^^'^ distribution. Since the chosen M^j in- 
tervals significantly exceed the mass resolution for all 
masses, we do not correct the mass spectrum for reso- 
lution effects. 



The measured cross section for e+e^ — + AA is shown in 
Fig. 13 and listed in Table IV. The quoted errors are sta- 
tistical and systematic. The latter includes the system- 
atic uncertainty in detection efficiency, the uncertainty 
in total integrated luminosity (1%), and the uncertainty 
in the radiative correction (1%). The only previous mea- 
surement of the e+e^ AA cross section, IOOI35 
2.386 GeV [6], is in agreement with our results. 



The e+e~ — ^ ylyl cross section is a function of two form 
factors. Due to the poorly determined \Ge/Gm\ ratio 
they cannot be extracted from the data simultaneously 
with reasonable accuracy. We introduce an effective form 
factor (Eq.(4)) which is a linear combination of jG^p and 
|GmP- The calculated effective form factor is shown in 
Fig. 14 and hsted in Table IV. 
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FIG. 12: The mass dependence of the FIG. 13: The e+e" -> AA cross section FIG. 14: The measured A effective form 

mass resolution calculated as the RMS measured in the present experiment com- factor, 

of the M^x - ^ax" distribution in MC pared to the DM2 [6] measurement, 
simulation. 



F. J/ip and rp{2S) decays into A A 

The differential cross section for ISR production of a 
narrow resonance (vector meson V), such as J/ip, decay- 
ing into the final state / can be calculated using [20] 



dais, B*) _ l2n^T{V e+e-)B{V /) 



d cos 9* 



my s 



(11) 

where my and r(V^ — > e~^e ) are the mass and electronic 
width of the vector meson V, xq = 1 — my/s, and B{V 
/) is the branching fraction of V into the final state /. 
Therefore, the measurement of the number of J/ip — s- AA 
decays in e+e~ — s- AAj determines the product of the 
electronic width and the branching fraction: T{J/tJj — > 
e+e-)B{J/ip^ AA). 

The AA mass spectra for selected events in the J/ ifj and 
ip{2S) mass regions are shown in Fig. 15. We determine 
the number of resonance events by counting the events 
in the signal region indicated in Fig. 15, and subtracting 
the number in the two sidebands. The following numbers 
of J/ip and ip{2S) events arc obtained: iVj/^ = 142 ± 12 



and N, 



17 ±4. A possible background due to ip 



ppTT~^'K~ decay is estimated using the two-dimensional 
distribution of the masses of A and A candidates. It 
is found to be 0.5lQ'g events for J ftp and negligible for 
^(25)- 

The detection efficiency is estimated from MC simu- 
lation. The event generator uses the experimental data 
for the angular distribution of the A in J/tp ^ AA de- 
cay. This distribution is described by 1 4- a cos^ 6 with 
a = 0.65 ± 0.010 [14-16]. For the ip{2S) the value 
a ~ 0.69 predicted in [21] is used. The error in the 
detection efficiency due to the uncertainty of a is negli- 
gible for the J/ip and is taken to be 5% for the '0(25'). 
The efficiencies corrected for data-MC simulation differ- 
ences are 0.062±0.004 for the J/ip and 0.059±0.005 for 
the ip{2S). 



The cross section for e~^e 
y < 160° is calculated as 

cr(20° <e*< 160°) 



iPj AA-f for 20° < 



eRL' 



(12) 



yielding (9.8 ± 0.9 ± 0.6) fb and (1.2 ± 0.3 ± 0.1) fb for 
the J/ip and the ip{2S), respectively. The radiative- 
correction factor R = a/aBom is 1.007 ± 0.010 for the 
J/ip and 1.011 ± 0.010 for the ip{2S), obtained from MC 
simulation at the generator level. 

The total integrated luminosity for the data sample is 
(230 ± 2) fb~^. From the measured cross sections and 
Eq. (11), the following products are determined: 



r(j/^ 
r(0(25) - 



^ e+e-)B{J/ip - 
e+e~)B{ip{2S) 



AA) = (10.7 ±0.9 ±0.7) eV, 
AA) = (1.5 ±0.4 ±0.1) eV. 



The systematic errors include the uncertainties in detec- 
tion efficiency, integrated luminosity, and the radiative 
correction. 

Using the world-average values of the electronic 
widths [13], the ip AA branching fractions are cal- 
culated to be 



B{J/ip AA) = (1.92 ± 0.21) X 10-^ 
B{ip{2S) AA) = (6.0 ± 1.5) x 10""* 



Both results are higher than the current world-average 
values [13]: (1.54 ±0.19) x 10"^ ^nd (2.5 ± 0.7) x lO"'*, 
but in reasonable agreement with the more precise recent 
measurements: (2.03 ± 0.15) x 10"^ by BES [16] and 
(3.33 ± 0.25) X 10-4 by CLEO [22] and BES [23]. 



G. Measurement of the A polarization 

A nonzero relative phase (p between the electric and 
magnetic form factors leads to polarization of the outgo- 
ing baryons. The exact formula for the A (A) polarization 
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FIG. 15: The AA mass spectra in the mass regions near the J/^ (a) and the '4>{2S) (b). The arrows indicate the boundaries 
between the signal regions and sidebands. 




FIG. 16: The distribution of C,fmax for FIG. 17: The distribution of cos Op^ for selected e^e AAy events with j < 
selected simulated events for e+e^ 2.8 GeV/c^ in simulation (a) and in data (b). 

AA'y with M^j < 2.8 GeV/c^ 



vector (^f is given in the Appendix. The polarization is 
proportional to sin 0. The magnitude of the polarization 
Cfmax = Cfi4' = ■"'/S) calculated under the assumption 
that \Ge\ = |Ga/| for simulated e+e^ AA'y events 
with M^j < 2.8 GeV/c^ is shown in Fig. 16. The sim- 
ulated events were reweighted according to the AA mass 
spectrum observed in data. The average value of Cfmax 
is equal to 0.285. The A polarization can be measured 
using the correlation between the direction of the A po- 
larization vector and the direction of the proton from A 
decay: 

diV" 

— — = A(l + a^C/cos0K), (13) 

U cos np(^ 

where is the angle between the polarization axis 
and the proton momentum in the A rest frame, and 
aA = 0.642 ± 0.013 [13]. For A, aj = -a^. The dis- 
tribution of cos0p^ for simulated e+e^ AAj events 
with M^j < 2.8 GeV/c^ (there is no A polarization in 
the simulation) is shown in Fig. 17(a). We combine the 



A and A distributions taking into account the different 
signs of a A and aj. Since the distribution is flat, we 
conclude that there is no dependence of the detection ef- 
ficiency on cos 0pc^ . A fit to the distribution using a linear 
function gives slope consistent with zero. 

The same distribution for data is shown in Fig. 17(b). 
In each angular interval the background is subtracted us- 
ing the procedure described in Sec. IIIB. The distribu- 
tion is fitted using a linear function. The slope is found 
to be 0.020 ± 0.097. The corresponding symmetric 90% 
CL interval for A polarization averaged over the AA mass 
range from threshold to 2.8 GeV/c^ is 

-0.22 < C/ < 0.28. 

Under the assumption \Ge\ = \Gm\ {Cfmax = 0.285), 
which does not contradict the data, this interval can be 
converted to an interval for sin^ as follows: 

-0.76 < sin0 < 0.98. 
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Our statistics allow only very weak limits to be set on 
sin (j). 



IV. THE REACTION e+e" S° S^-y 
A. Event selection 

The Z"" hyperons are detected via the decay Z"" A'^ 
(the branching fraction is 100% [13]). Therefore, the pre- 
liminary selection of e+e^ U'^S^j candidate events is 
similar to that for e^e^ — > AA-y. In addition, we re- 
quire that an event contain at least two extra photons 
with energy greater than 30 MeV. To suppress combina- 
torial background from events not containing two A^s in 
the final state, we apply a tighter selection criterion on 
the mass of a yl (A) candidate: 1.110 < Af^^- < 1.122 
GeV/c2. 

For events passing the preliminary selection, we per- 
form a kinematic fit to the e^e~ AA'j^fj hypothe- 
sis. The photon with highest E* is assumed to be the 
ISR photon. The fitted momenta of two other photons 
and yl-baryons are used to calculate A'y and A'y invariant 
masses. For S'^ and candidates these masses must 
to be in the range 1.155-1.23 GeV/c^ (the nominal 
mass is 1.192642(24) GeV/c^ [13]). We require that an 
event contain at least one E° and one 17° candidate. For 
events with more than three photons we iterate over all 
possible photon combinations and find the one contain- 
ing S'^ and Z'" candidates and giving the lowest for 
the kinematic fit. 

The distribution of the of the kinematic fit (x'ss) 
for simulated e+e" S'^S'-'j events is shown in Fig. 18. 
We select data events with Xss < 20 for further analysis; 
as before, a x^ control region (20 < Xss < ^0) is used 
for background estimation and subtraction. 

To suppress background from e+e^ AA'y and 
e+e^ AS^'j events with extra photons we also perform 
kinematic fits to the ylyl7 and AS'^j hypotheses. The 
AS°j fit is a fit to the e+e" AA-y-/ hypothesis. The 
photon with highest E* is assumed to be the ISR photon. 
The other photon taken in combination with the A or yl, 
must give an invariant mass value in the range 1.155- 
1.23 GeV/c^, where the mass is calculated using fitted 
momenta. The x'aa distributions for simulated events 
corresponding to e'^e" AA'y and e+e^ S'^S'^j are 
shown in Fig. 19. The requirement x'aa > 20 rejects 93% 
of AA'y events and only 3% of signal events. Similarly, the 
Xas distributions for simulated events for e+e~ — > AS'^'y 
and e~^e~ —* S^S^j are shown in Fig. 20. The require- 
ment Xas > 20 again rejects 93% of AU^^ events, but in 
this case removes 30% of the signal events. Data events 
with x'as < 20 are used to estimate the level of AS'^^ 
background. 

The scatter plots of the invariant mass of the 17° 
candidate versus the invariant mass of the S'^ candi- 
date for the selected data events and simulated e'*"e~ 



r°r°7 events are shown in Figs. 21(a) and (b), respec- 
tively. Of the two possible A'y and A'y combinations, 
we plot only the combination with the smaller value of 
{Maj — + {^A-t ~ "^^)^' where ms is the nom- 

inal 17° mass. The 17° 17° mass spectrum for the data 
events with the additional requirement that the 17° and 
17° candidate mass values satisfy 1.180 < Ma^ < 1.205 
GeV/c^ (central box in Fig. 21(a)), is shown in Fig. 22. 
An excess of signal events is seen at masses below 3.0 
GeV/c^. There are also about 20 events near the J ftp 
mass, corresponding to J/ip 17° 17° decay. The two 
events near 3.7 GeV/c^ may be due to ^(23) ^ 17° r° 
decay. The mass distribution for 17° and 17°candidates 
from the J ftp region is shown in Fig. 23. The spectrum 
is obtained as the difference of the spectrum from the 
region 3.05 < M^^ < 3.15 GeV/c^ and that from the 
sideband region (3.00-3.05 and 3.15-3.20 GeV/c^). We 
see that simulation reproduces the 17° lineshape quite 
well. 



B. Background subtraction 

Background processes can be divided into three classes, 
namely those with zero (e+e" AA'y, AAtt'^, AAtt'^j, 
... ), one {e+e" AS°-/, AS^n°, yir°7r°7, . . . ), and 
two 17°'s (e+e- -> 17°17°7r°, 17°17°7r°7, ... ) in the final 
state. To separate events with two 17° 's from events with 
no 17° 's and one 17° we use the differences in their two- 
dimensional distributions of invariant mass values of the 
17° and 17° candidates. 

Background events from e+e^ 17°17°7r° with an un- 
detected low-energy photon or with merged photons from 
7r° decay yield a low value of x^ when reconstructed under 
the 17°17°7 hypothesis, and can not be separated from 
the process under study. Special selection procedures are 
applied in order to estimate this background. The proce- 
dures are similar to those used to study background from 
e+e- AAtt° in Sec.IIIB. No 17°r°7r° candidates are 
found in data, and we estimate that the background from 
this process does not exceed 5 events. Assuming that the 
dibaryon mass spectrum in e+e~ 17°17°7r° is similar to 
that for e+e~ pp7r° [1] we find that about 70% of the 
e+e~ 17°17°7r° events are expected to have 17°17° mass 
less than 3 GeV/c^. Two- 17° background other than that 
from 17°17°7r° can be estimated using the difference in the 
distributions for signal and background events. 

The 3x3 two-dimensional histograms of My^^ vs A/j^ 
(dashed lines in Fig. 21) for events from three classes: 

Xls < 20, xIe > 20, 
20 < xIe < 40, Xae > 20, 
xIe < 20, xls < 20, 

are fitted simultaneously. The second histogram is used 
to determine two- 17° background. From the third his- 
togram we estimate the e+e~ — s- yll7°7 background. 
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FIG. 18: The Xee distributions for sim- 
ulated events for e'''e~ — > S'^E'^'y. 



FIG. 19: The Xaa distributions for simu- 
lated events for e'^e~ AA'y (solid his- 
togram) and e^e~ — > S^E^^ (dashed 
histogram). 



FIG. 20: The Xj\.s distributions for sim- 
ulated events for e^e~ — > AS^'y (solid 
histogram) and e^e~ E^E^j (dashed 
histogram) . 
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FIG. 21: Scatter plots of the invariant mass of the X'" candidate versus the invariant mass of the X''' candidate for selected 
data events (a), and simulated e+e~ — > E" E'^j events (b). 
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FIG. 22: The X'^X'" mass spectrum for the FIG. 23: The A"f and Aj invariant mass 

selected data events. distribution for data (points with error 

bars) and simulated (histogram) events 
from the J/tp region. 
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Each histogram is fitted using the following function: 



N2flf + Niflf + Noflf, 



'1,3 



where A'o, iVi, and N2 are the numbers of events with 
zero, one, and two 17" 's in the final state. The func- 
tions f^^ and /^^ are taken from e+e" AE'^j and 
e+e^ Z'°Z'°7 simulations. The probability density 
function for zero- events is the product of two identical 
linear functions of Maj and Mj^ and a function taking 
into account the correlation between masses of the 17° 
and Z!^ candidates. This last function is extracted from 
e~^e~ — *■ AA'y simulation. The correlation arises from 
our choice of one of the two possible combinations of A 
and A candidates with photons, and is about 15% for the 
central mass bin. 

In order to find the number of signal events and esti- 
mate the background we use the following relations: 

N^ixls < 20, Xae > 20) = N2s + N^b, 
iV2(20 < xIe < 40, Xas > 20) = a^N^s + fiiN^i,. 
N2{x%E < 20, Xae < 20) - a2N2s + P2N2b. 

Niixls < 20, Xas > 20) = iVi. + Nib, 
iVi(20 < xIe < 40, xls > 20) = 71 ^is + ^lA^ 



lb- 



7Vi(xl;^<20, s < 20) = 72iVi. 



-S2N 



lb, 



where N2S is the number of signal S'^S'^j events and 
N2b is the number of two-Z"*^ background events in the 
signal region {xss < 20, Xas > 20), Nis is the num- 
ber of ylZ'°7 events and Nib is the number of one- 17° 
events from all other processes in the signal region; N2S, 
N2b, Nis, and Nib, are then free fit parameters. The 
coeflScients a^, (3i, 7;, and Si are obtained from simula- 



tion of e+e" A.'yj^'j, 

AE^j, e+e~ AS^tt^j, respectively. For the coeffi- 
cients most critical to the analysis, ai = 0.22 ± 0.03 and 
/?! = 1.5 ±0.3, the errors include uncertainties due to the 
data-MC difference in the x^ distributions for the kine- 
matic fits. The other 6 free parameters are the numbers 
of zero-Z° events in the three histograms, and the slopes 
of the linear functions describing the mass distributions 
for these events. 

The fit results for Z'°i7° masses below 3 GeV/c^ are 
shown in Table V, together with the predictions from 
JETSET simulation. The one-S° background is domi- 
nated by the e+e" AS^j process. The number of 
one- 17° events from other processes is found to be con- 
sistent with zero. The numbers of e'^^e^ — > S^S'^j and 
e+e" AS°"f events with Xae < 20 are 7.7^11 and 
15.3+7;y, respectively. 

The fitting procedure was performed in five S'^S'^ mass 
ranges, and the number of signal events found in each is 
listed in Table VII. For S°S° masses below 3 GeV/c^ 
we observe an excess of signal events over background. 
The significance of the observation of E^E^ production 
in the mass region below 3.0 GeV/c^ is 2.9(t. For other 
mass bins we list upper limits at the 90% CL. 
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FIG. 24: The X'''X'° mass spectrum for the mass region near 
the J/ip. The arrows indicate the boundaries between the 
signal region and sideband regions. 



C. Cross section and form factor 

The cross section for e+e" E^E'^ is calculated from 
the E^E^ mass spectrum according to Eqs.(9-10). 

The detection efficiency is determined from MC simu- 
lation and then corrected for data-MC simulation differ- 
ences in detector response. The model dependence of the 
detection efficiency due to the unknown \Ge/Gm\ ratio 
is estimated to be 5% (see Sec. HID). The efficiency 
corrections summarized in Table VI were discussed in 
Sec. HID. On the basis of our analysis of ISR processes 
with photons in the final state [24], we enlarge the sys- 
tematic error in the correction for the x^ selection inter- 
val. The correction for trigger inefficiency is removed, 
since for e+e~ E'^E'^j, trigger inefficiency is less than 
0.001 both in data and in MC simulation. The corrected 
detection efficiencies are listed in Table VII. The over- 
all uncertainty in efficiency takes into account simula- 
tion statistical error, model uncertainty, the error in the 
A piT^ branching fraction, and the uncertainty of the 
efficiency correction. 



The measured values of the e+e 



E°E° 



cross sec- 
tion are listed in Table VII, together with the effective 
form factor values calculated according to Eq.(4). The 
quoted errors on the cross section are statistical and sys- 
tematic. The latter includes systematic uncertainty in 
detection efficiency, the uncertainty in total integrated 
luminosity (1%), and radiative correction uncertainty 
(1%). This is the first measurement of the e+e~ — > E'^E'^ 
cross section. The upper limit set by DM2 [6] at 2.386 
GeV ( < 120 pb) is consistent with our measurements. 



D. J/V^ decay into 

The E^E'^ mass spectrum for selected events in the 
J/i/j mass region is shown in Fig. 24. We determine the 
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TABLE V: Comparison of the fit results for 17° 17° masses below 3 GeV/(? and tlie predictions from JETSET simula- 
tion; A'2s, No, Ni, N2b are the fitted numbers of signal, zero-, one-, and two- 17" background events in the signal region 
iXss < 20, Xas > 20), respectively, A'^^o ^o^o is expected number of background events from e'*'e~ — > X'^X'^tt" process. 





N2s 


No 


Ni 


N2b 




data 




11.3 ±4.8 


1.2 ±0.6 


2.8 ±4.8 


< 4 


JETSET 


33 ±5 


3.1 ± 1.4 


1.2 ±0.8 


< 1.4 


0.6 ±0.6 



TABLE VI: The values of the various efficiency corrections 



for the process e e 



17° r^. 



effect 


s., (%) 


xijE < 20 


-2.0 ±6.0 


track reconstruction 


-1.0 ±3.8 


p nuclear interaction 


+1.0 ±0.4 


PID 


+0.6 ±0.6 


photon inefficiency 


-3.9 + 0.9 


photon conversion 


+1.2 + 0.6 


total 


-4.1 + 7.2 



800 - 




400 - 



number of resonance events by counting the events in 
the signal region indicated in Fig. 24 and subtracting the 
number in the two sidebands. The net number of J/ip 
decay events is then 30 ± 6. 

The detection efficiency is estimated from MC simu- 
lation. The event generator uses the experimental data 
on the angular distribution in J/i/; — > S'^E^ decay, 
whichisl + Q;cos2t/ieiawitha= -0.1±0.2 [14-16]. The 
error in the detection efficiency due to the uncertainty in 
a is negligible. The efficiency corrected for data-MC sim- 
ulation differences is e = 0.022±0.002. 

Using Eqs. (11,12) the following product is determined: 

TiJ/tP e+e-)BiJ/i' = (6.4 ± 1.2 ± 0.6) eV. 

The systematic error includes the uncertainties in detec- 
tion efficiency, integrated luminosity, and in the radia- 
tive correction. Using the PDG value of the electronic 
width [13], the J/i/' — > i7°Z'° branching fraction is calcu- 
lated to be 

B{J/ip ^ (1.16 ± 0.26) X IQ-^ 

Our result is in agreement with the world average value 
(1.31 ±0.10) X 10-3 [13]. 

We also observe 2 events in the •ip{2S) region with 
zero background, estimated from the sidebands. This 
number agrees with the 2.5 ± 0.4 events expected from 
the measured branching fraction B{ip{2S) Z'"Z'°) = 
(2.51 ±0.31) X 10-4 [22, 23]. 



FIG. 25: The Xae distributions for simulated e+e- AE^-y 
events. 



V. THE REACTION e+e" AS"-/ 
A. Event selection 

The preliminary selection of e'^e~ — > AS'^j events is 
similar to that for e~^e~ — + ylyl7. Additionally we require 
that an event candidate contain at least one extra photon 
with energy greater than 30 MeV. To suppress combina- 
torial background from events not containing two /I's in 
the final state we require that the mass of the A (A) 
satisfy 1.110 < Mp^- < 1.122 GeV/c^. 

For events passing the preliminary selection, we per- 
form a kinematic fit to the e+e- — > AS'^j hypothesis as 
described in Sec. IV A. The Xas distribution for sim- 
ulated AS'^j events is shown in Fig. 25. We select the 
events with Xas < 20 for further analysis. The control 
region (20 < Xas < 40) is used for background estima- 
tion and subtraction. To suppress background result- 
ing from e~'"e- —f AAj events with an additional photon 
we perform a kinematic fit to the ylyl7 hypothesis and 
require that Xaa > 20. The Xaa distributions for simu- 
lated e+e" — > AA'y and e+e" AS'^j events are shown 
in Fig. 26. The Xaa > 20 cut rejects 95% of the AA-f 
events at the cost of 20% of the signal events. 

The distribution of 17° candidate invariant mass for 
data events passing the AE'^j selection process is shown 
in Fig. 27. For each event we plot only the A{A)j com- 
bination closer to the nominal i7° mass. The AE'^ mass 
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TABLE VII: The S°E° invariant mass interval (M^jj), net number of signal events {Ns 



luminosity (L), measured cross section (a), and effective form factor {F) for e"* 
statistical and systematic. For the form factor, the total error is listed. 



detection efflciency (e), ISR 
The quoted errors on a are 



(GeV/c2) 




e 


L 

(pb-i) 


a 
(pb) 


\F\ 


2.385-2.600 
2.600-2.800 
2.800-3.000 
3.200-3.600 
3.800-5.000 


10.3+1-^ 

lA+_li 

< 2.3 

< 2.3 


0.024 ± 0.002 
0.025 ± 0.003 
0.026 ±0.003 
0.023 ±0.003 
0.023 ±0.002 


14.3 
14.7 
16.1 
39.9 
180.4 


30 ± 13 ± 3 

3.4tf | ± 0.4 

< 2.5 

< 0.5 


047+°"^" 

"■'J^' -0.017 

n 021+°°^* 

U.UZi_Q 

< 0.019 

< 0.011 





M^^(GeV/c^) 



FIG. 26: The Xaa distributions for simulated e"'"e~ — > AA'y 
events (solid histogram) and e+e^ — > AS^^ events(dashed 
histogram) . 



distribution for selected data events with invariant mass 
of the r° candidate in the 1.185-1.205 GeV/c^ range is 
shown in Fig. 28. We expect that the e+e" S°IJ°-f 
process results in a significant contribution to the selected 
event sample. In particular, the peak in the AS'^ mass 
spectrum near 3 GcV/c^ is from J/tp ^ i7°i7° events 
with a missing or excluded photon. 



B. Background subtraction 

The background processes can be divided into three 
classes, namely those with zero (e+e^ ^^7; AAtt^, 
AAn^-f, ... ), one {e+e' -> AS^tt^, AS°n°-f, ... ), and 
two IJ°'s (e+e- S^S^-f, S^S^tt^, S^S^Ti^-f, ... ) in 
the final state. To separate onc-S^ events from events 
with no E'^ we use the difference in mass distribution for 
the respective (^°) candidates. 

To determine two-S^ background we select a clean 
sample of two-Z"*^ events. To do this the E'^S'^j crite- 
ria (Sec. IV A) are used with the additional requirements 
that 1.180 < Maj < 1.205 GeV/c^ for the E° and E° 
candidates, and that Xas < 100. The latter requirement 



FIG. 27: The distribution of the invariant mass of the S and 
X'" candidates for the selected AS'^j candidate. 
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FIG. 28: The AE'^ invariant mass spectrum for data events 
with the invariant mass of the S° candidate in the 1.185-1.205 
GeV/c^ range. 



is needed to obtain a useful distribution. Using 

the ratio of detection efficiencies for two-Z"*^ and AE^j 
selections (k = 0.80 ± 0.05), we can convert the num- 
ber of events in the two- 17° sample to an estimate of the 
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number of background events in the AE^j sample. 

The background from e~^e~ — > AE'^tt'^ events with an 
undetected low-energy photon or with merged photons 
from TT*^ decay cannot be separated from the process un- 
der study. The experimental data with special selection 
criteria are used to estimate this background. The proce- 
dure is similar to that used in the study of e+e" — > AAtt^ 
background in Sec. Ill B. We selected two E^E'^tt^ can- 
didates with an expected background from e~^e~ — > AAj 
and e+e" AE^j processes of 0.5 ±0.2 events. To sup- 
press the AE'^tt'^ background in the AE'^^ sample we re- 
ject events with 0.10 < < 0.17 GeV/c^, where 
is the invariant mass of the most energetic photon in an 
event and another photon with energy greater than 0.1 
GeV. This removes about 1/3 of AE^tt'^ events and less 
than 1% of signal events. After applying this selection 
criterion the rates at which e+e^ — > AE'^tt^ events are 
selected a&_yir°7 or AE°n° are in the ratio (2.1 ± 0.2), 
and the AE'^tt'^ background in the AE'^j event sample is 
estimated to be (3.1 ± 2.2) events. We assume that the 
dibaryon mass distribution for the e~^e~ — > AE'^tt'^ pro- 
cess is similar to that for the e'^e^ ppir^ process [1]. In 
particular, about 70% of the e+e" AE'^tt^ events have 
AE° mass less than 2.9 GeV/c^. Both observed AE^n^ 
events lie in this mass region. 

The one-i7'' backgroimd other than AE^tt^ is esti- 
mated using the difference in the distributions for sig- 
nal and background events. Two histograms of Ma^ for 
events with x\x: ^ ^0 and with 20 < Xae < ^0 arc fitted 
simultaneously to the sum of the distributions for signal 
and background 



n, = NlHu+N2H2^ + NoH^ 



Oi, 



(14) 



where A^i, iV2, and iVo arc the numbers of events con- 
taining one, two, and zero E^'s in the final state, re- 
spectively. The onc-E^ events are the signal events 
with a possible contribution from the background pro- 
cesses AE'^n^, AE'^n^j, etc. The function Hi describ- 
ing the mass distribution of one-Z''^ events is calculated 
using e+e^ —^ AE'^j simulation. The distribution of 
two-E*^ events is taken from e+e~ — > E'^E^j simula- 
tion. The parameter iV2 is fixed by addition of the term 
— ln/p(no; fJ-o) to minimize the likelihood function. Here 
fp is a Poisson distribution, uq is number of events in 
the two- 17" sample described above and /io = KrN2. The 
scale factor n = 0.80±0.05 is found from e+er E°E°-/ 
simulation as the ratio of detection efficiencies for two- 17" 
and AE'^'j selections. The factor r « 1.1 takes into ac- 
count the purity of the two-Z"" sample, which is (90±5)%. 
It should be noted that the two-Z"" sample contains not 
only E'^E^j but also events from the processes e+e~ 
Z^Z^tt", E'^E'^tt^j, etc. A similar approach is used to in- 
troduce the AE'^tt'^ background into the fit. The shape of 
the zero-i7° background (Hq) is modeled using the mass 
distribution for the e~'"e~ —> AAj process. The distribu- 
tion is parametrized as Hq = (l-|-ao(Myi^ — mi;))/(M^), 
where m^; is nominal E'^ mass. The function f(Mj[f) 
describes the deviation from a linear function due to our 
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FIG. 29: The distribution of the invariant mass of the 1!° {1!°) 
candidate for data events with M^j^ < 2.9 GeV/c^ (points 
with error bars). The solid histogram shows the result of the 
fit described in the text. The dotted histogram shows the 
contribution of zero- 17° background. The difference between 
dashed and dotted histograms is the contribution of two-X"" 
background. 
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FIG. 30: The distribution of selected data events (points with 
error bars) over chosen mass intervals. The histogram shows 
fitted background. 



choice of one of the two A{A)j combinations. This func- 
tion is equal to unity at the end points of the mass in- 
terval, and is about 2 at the center. We checked that 
the function Hg with oq as free parameter provides a 
good description of the mass distributions for simulated 
e+e^ AAtt'^ and e+e^ AAtt'^j events, and data ylyl7 
events selected by requiring Xaa < ^0 and Xas < 20. 

The one- 17° background other than that from AE^tt^ 
is estimated from the fit according to Eq.(6). The coef- 
ficients Psig and Pbkg arc obtained from the signal and 
ylZ'°7r°7 simulation and take the values 0.15 ± 0.02 and 
1.5 ± 0.3, respectively. Their errors are enlarged to take 
into account the data-MC simulation difference in the 
distributions resulting from the kinematic fits. 
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The fit results for ylZ'" masses below 2.9 GeV/c^ are 
shown in Fig. 29 and summarized in Table VIII, together 
with the predictions from the JETSET simulation. 

The fitting procedure was performed in eight AS^ mass 
intervals, and the resulting data distribution is compared 
to the fitted background in Fig. 30. An excess of signal 
events over background is seen only for AS^ masses below 
2.9 GeV/c^. The number of signal events in each mass 
interval is listed in Table X; 90% CL upper limits are 
given for the intervals with M^^ > 2.9 GeV/c^. The 
significance of the observation of AS'^ production in the 
mass region below 2.9 GeV/c^ is 3.3(t. 



C. Cross section and form factor 

The cross section for e+e~ AS'^ is calculated from 
the AS'^ mass spectrum according to Eqs.(9-10). 

The detection efficiency is determined from MC simu- 
lation and then corrected for data-MC simulation differ- 
ences in detector response. The model dependence of the 
detection efficiency due to the unknown \Ge/Gm\ ratio 
is estimated to be 5%. The efficiency corrections summa- 
rized in Table IX were discussed in Sees. HID and IV C. 
The corrected detection efficiencies are listed in Table X. 
The uncertainty in efficiency takes into account simula- 
tion statistical error, model uncertainty, the error on the 
A — > p7r~ branching fraction, and the uncertainty in the 
efficiency correction. 

The measured values of the e~^e^ AS'^ cross section 
are listed in Table X, together with those of the effec- 
tive form factor. The quoted cross section errors are 
statistical and systematic. The latter includes system- 
atic uncertainty in detection efficiency, the error on the 
total integrated luminosity (1%), and the radiative cor- 
rection uncertainty (!%)• This is the first measurement 
of the e+e" — > AE'^ cross section. The upper limit set 
by DM2 [6] at 2.386 GeV ( < 75 pb) is consistent with 
our results. 

Assuming that all events in the 2.90-3.30 GeV/c^ mass 
range result from J/tp ^ AI!° decay we obtain an upper 
limit for the J/ip ^ AS'^ branching fraction B{J/ip 
AS") < 2 X 10"'*, which is slightly higher than the only 
other estimate, S(J/'(/' ^ AS") < 1.5 x lO""* [26]. 



VI. SUMMARY 

The processes e+e~ AAj, AS"j, and E"S"j 
have been studied for dibaryon invariant mass up to 5 
GeV/c^. From the measured dibaryon mass spectra we 



^ For the e+e~ — > AE'^ process, Eq. (3) must be modified by 
the substitutions /3 = (1 — {rriA — ms)^ /m^)2P^/m and r = 
m^/(mj\ + iTijj)^ [25], where is the baryon momentum. 
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FIG. 31: The measured dependence of the baryon form factors 
on dibaryon invariant mass. The proton data are taken from 
Ref. [1]. 

obtained the e+e~ AA, AS", and E"E" cross sec- 
tions and baryon effective form factors. Our results on 
the measurements of the various baryon form factors for 
dibaryon invariant masses above AA threshold are shown 
in Fig. 31. 

For e+e~ AA^y wc analyzed the A angular distri- 
butions in the mass range from threshold to 2.8 GeV/c^ 
and extracted the \Ge/Gm\ ratio. Our results are 

\Ge/Gm\ = 1.73l°J? for 2.23-2.40 GeV/c^, 

\Ge/Gm\ = 0.71+°;?? for 2.40-2.80 GeV/c^, 

and are consistent both with \Ge/Gm \ = 1 and with the 
results for e+e" — > pp [1], where this ratio was found to 
be significantly greater than unity near threshold. 

The measurement of the A polarization enables the 
extraction of the relative phase between the A electric 
and magnetic form factors. The limited statistics of the 
present experiment allow us to set only very weak limits 
on this phase: 

-0.76 < sin0 < 0.98. 

From the events in the J/V' and ipi^S) regions the 
products, 

T{J/xP e+e-)B{JMj AA) = (10.7 ± 0.9 ± 0.7) cV, 
r(J/i/) ^ e+e-)B{J/Tp ^ E"E") ^ (6.4 ± 1.2 ± 0.6) eV, 
F(V'(25') ^ e+e-)B{ilj{2S) ^ AA) = (1.5 ± 0.4 ± 0.1) cV, 

have been measured, and, using the known e+e~ partial 
widths, the corresponding branching ratios have been ob- 
tained: 

B{J/%p -> AA) = (1.92 ± 0.21) X 10"^ 
B{J/iIj E"e") = (1.16 ± 0.26) x 10-^ 
^(^(25*) ^AI) = (6.0 ± 1.5) X 10"''. 



23 



TABLE VIII: Comparison of the fit results for AS^'' masses below 2.9 GeV/c^ and the predictions from JETSET simulation; 
Nis, No, Nib, N2 and Ny^^o^o are the fitted numbers of signal, zero-, one-, two-X"", and X'^X'^'tt" background events with 
Xas < 20, respectively. 





Nis 


N2 


No 






data 


24.1 ±8.4 


13.8 ±4.4 


17.0 ±7.8 


<5 


3.1 ±2.2 


JETSET 


50 ±6 


17 ±3 


3.0 ± 1.5 


0.6 ±0.6 


1.2 ±0.9 



TABLE IX: The values of the various efficiency corrections 
for the process e"'"e~ AS'-^'j. 



effect 


s,, (%) 


Xae < 20 


-2.0 ±6.0 


track reconstruction 


-1.0 ±3.8 


p nuclear interaction 


±1.0 ±0.4 


PID 


±0.6 ±0.6 


photon inefficiency 


-2.6 ±0.6 


photon conversion 


±0.8 ±0.4 


total 


-3.2 ±7.2 
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APPENDIX: ANGULAR DISTRIBUTIONS AND A POLARIZATION IN THE e+e" AA-y REACTION 

The formulae given in this section arc taken from Rcf. [3]. The process e~^e~ AAj is considered in the e^e^ 
center-of-mass frame, where the electron has momentum p and energy e, and the photon has momentum k and energy 
uj. The A momentum P is given in the AA rest frame. The differential cross section summed over the polarization of 
one of the final particles is given by 

= W.^Sq^il'^il-ur] ^^' + Cfs), 2\Gm\M + N^) + {^\Ge\' ~ IGmP) ([n x f]^ + [N x f]^). 



Cf = ^Im(Gl,GM) ((n • f)[n x f] ± (N • f)[N x f]) ; n=^,u^^,N= - + - ^)(- ^ 

1 2£-CJ 



iV2 = (n.i.)2±-^, 7=^_, Q^^fli^), P=|P| = ^QV4-m2, f = -. 
Here s and C,f arc the spin and polarization vectors of the A in its rest frame. 
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